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Abstract 

We develop a new method for simultaneous timing and spectral studies of Gamma Ray Burst 
(GRB) prompt emission and apply it to make a pulse-wise description of the prompt emission 
of GRB 090618; the brightest GRB detected in the Fermi era. We exploit the large area (and 
sensitivity) of Swift /BAT and the wide band width of Fermi /GBM to derive the parameters for 
a complete spectral and timing description of the individual pulses of this GRB, based on the 
various empirical relations suggested in the literature. We demonstrate that this empirical model 
correctly describes the other observed properties of the burst like the variation of the lag with 
energy and the pulse width with energy. The measurements also show an indication of an increase 
in pulse width as a function of energy at low energies for some of the pulses, which is naturally 
explained as an ofF-shoot of some particular combination of the model parameters. We argue 
that these model parameters, particularly the peak energy at the beginning of the pulse, are the 
natural choices to be used for correlation with luminosity. The implications of these results for 
the use of GRBs as standard candles are briefly described. 



Subject headings: Gamma-ray burst: general - 
data analysis — Methods: observational 

1. INTRODUCTION 



- Gamma-ray burst: individual (GRB 090618) — Methods: 



The phenomenon of Gamma Ray Burst (GRB) 
has raised many unresolved issues. For example, 
the nature of the compact object in the progeni- 
tor responsible for the jet emission is not yet un- 
derstood. The widely accepted scenario for the 
GRB emission, supported by strong observational 
evidences (van Paradijs et al. 2000), is that of a 
catastrophic energy release from a highly massive, 
rapidly rotating, low metallicity star towards the 
end of its life. This catastrophic process for the 
long GRB class (Tgo 2 s, where Tgo is the time 
taken by the burst to accumulate between 5% to 
95% of its prompt emission - see Kouveliotou et 
al. 1993; Norris et al. 2005) is attributed to the 
coUapsar or hypernova scenario (Woosley 1993; 
Paczynski 1998; Fryer et al. 1999; Meszaros 2006) 
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while the short GRBs (Tgo < 2 s) are thought 
of as the outcome of merging neutron star bina- 
ries or neutron star-black hole binaries (Paczynski 
1986; Meszaros 1992;1997; Rosswog et al. 2003a; 
b). Though the collapsar model supports the for- 
mation of a stellar mass black hole, it is also sug- 
gested that the compact object could be a magne- 
tar (Metzger et al. 2007; Metzger 2010). 

The prompt emission of long GRBs shows a wide 
range of structures in their light curve, though 
the time integrated energy spectra are adequately 
fit by a four parameter Band model (Band et 
al. 1993). Attempts have been made to fit the 
light curve of the prompt emission with empirical 
models such as the Fast Rise Exponential Decay 
(FRED) model (Kocevski et al. 2003) and the ex- 
ponential model (Norris et al. 2005). These mod- 
els assume an underlying pulse structure of the 
prompt emission. The curve fitting shows that 
the pulses are nearly simultaneous in a wide range 
of energy bands (Norris et al. 2005) and have 
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self similar shapes (Nemiroff 2000). Hakkila et 
al. (2011) showed that the observed properties of 
the short pulses of long GRBs correlate the same 
way as typical pulses of short GRBs. All these 
properties of pulse emission indicate that despite 
being diverse in their many features, the prompt 
emission of GRB has a simple underlying mecha- 
nism of pulse emission which can be characterized 
by empirical models irrespective of the progenitor 
or the environment of its formation (Hakkila et al. 
2011). 

The properties of the prompt emission of GRBs 
like Epgak (the energy at which the emission peaks 
- Amati et al. 2002), the lag between the light 
curves at two energies (Norris et al. 2000), the 
average variability seen during the prompt emis- 
sion (Fenimore & Ramirez-Ruiz 2000), the rise 
time of the pulses (Schaefer 2002) are found to 
be correlated to the luminosity of the bursts and 
hence they are used as distance indicators (Schae- 
fer 2007; Gehrels et al. 2009). Most of these 
derived relations are empirical in nature and the 
analyses that are done to derive these parame- 
ters, however, are done separately in the time or 
the energy domain. For example, to derive i^peafe 
the time-integrated spectrum is used, ignoring the 
time variability and the pulse characteristics. Sim- 
ilarly, to derive lags, the time profile is exam- 
ined in a few broad energy bands. Attempts have 
been made to investigate the variations of these 
derived parameters in the other domain. For ex- 
ample, the variation of Ep^ak with time has been 
investigated by Kocevski & Liang (2003). A com- 
plete empirical description of GRB prompt emis- 
sion which describes the time as well as energy 
distribution, however, is lacking. This is probably 
due to the method of spectral parameter deter- 
mination which is forward in nature: that is one 
assumes a possible spectral model, convolves with 
the detector response and verifies the validity of 
the model by comparing with the data (Arnaud 
1996). Hence it is difficult to incorporate time 
variability parameters in this scheme. GRBs often 
have overlapping pulses which further complicates 
the extracting of model parameters for any unified 
empirical description. 

In this paper we develop a technique to give 
a complete pulse-wise empirical description of 
the prompt emission of a GRB using the time- 



integrated spectrum and the energy integrated 
light curve. We apply this method to the 
Fermi /GBM and the Swift /BAT data of GRB 
090618, which has the highest fluence among all 
438 GBM GRBs until March 2010 (Nava et al. 
2011). We exploit the large area of Swift/BAT to 
get accurate light curves and the large band-width 
of Fermi /GBM to get accurate energy spectra. 
We use the Band description for the time inte- 
grated spectral model (Band 2003) and the expo- 
nential model (Norris et al. 2005) for the energy 
integrated timing description for the pulses. We 
use the expression for the variation of Epeak with 
time given by Kocevski & Liang (2003) and de- 
velop a method to derive the model parameters 
(sec section 4.1.2). We demonstrate that these 
model parameters correctly describe the variation 
of pulse width with energy and the variation of 
the lag with energy. 

In §2 we give a brief description of GRB 090618, 

and in §3 we give the data analysis techniques 
including spectral analysis (§3.1), timing analy- 
sis (§3.2) and a pulse- wise joint {Swift /BAT and 
Fermi/GBM) spectral analysis (§3.3). In §4 we 
describe our method in detail, which gives a com- 
plete three dimensional pulse-wise description of 
the GRB. The model predicted timing results arc 
compared with the observations in §4.2 and the 
implications for global correlations are touched 
upon in §4.3. The results are discussed and the 
major conclusions are given in the last section (§5) . 

2. GRB 090618 

The brightest Gamma Ray Burst (GRB) in the 
Fermi era (as of March 2010), GRB 090618, was 
siimiltancously observed with Swift /BAT (Schady 
et al. 2009a; 2009b; 2009c), AGILE (Longo et 
al. 2009), Fermt/ GBM (McBreen et al. 2009), 
Suzaku/WAM (Kono et al.), Mnd/Konus-Wind, 
Coronas-Photon /Koims-KF (Golenetskii et al. 
2009), Coronas- Photon /RT-2 (Rao et al. 2009; 
2011). This is a long GRB (T90-II3 s) with a flu- 
ence of 3398.1 ± 62.0 x 10"'' erg/cm^ (integrated 
over At=182.27 s). The detection time reported 
by Swift /BAT is 2009 June 18 at 08:28:29.85 UT 
(Schady et al. 2009a) whereas the Fermi/GBM 
detection time is 08:28:26.66 UT (McBreen et al. 
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2009). This burst showed structures in the Ught 
curve with four peaks along with a precursor. A 
significant spectral evolution was also observed in 
the prompt emission which lasted for about 130 s 
from the trigger time (for a summary see Rao et 
al. 2011). 

The GRB afterglow was observed in the softer 
X-rays and optical regions. The optical after- 
glow was tracked by Palomar 60-inch telescope 
(Cenko 2009), Katzman Automatic Imaging Tele- 
scope (Perley et al. 2009) followed by various 
other optical, infrared, and radio observations. 
Observations using the Kast spectrograph on the 
3-m Shane telescope at Lick Observatory (Cenko 
et al. 2009) revealed a redshift of z = 0.54 (Cenko 
et al. 2009a). The X-ray afterglow was measured 
by the Swift /XRT (Schady et al. 2009b). Initially 
very bright in the X-rays, the flux rapidly decayed 
with a slope of ~6 before breaking after 310 s from 
the trigger time when the burst entered a shallower 
decay phase (slope 0.71±0.02 - Beardmore et al. 
2009). 

Rao et al. (2011) analyzed the prompt emis- 
sion of GRB 090618 using Coronas-Photon /RT-2 
and Swift/TiAT data. The light curves in vari- 
ous energy bands were fitted with a FRED pro- 
file to identify pulse characteristics such as width, 
peak position etc. A joint RT-2 and BAT spec- 
tral analysis of individual pulses was done us- 
ing the spectral model of Band. The reported 
time integrated low energy photon index (a), the 
high energy photon index (/3), and the peak en- 
ergy (Epeak) are -1.40, -2.50, and 164 keV, re- 
spectively. Ghirlanda et al. (2010) performed a 
time resolved spectral analysis of Fermi GRBs 
of known redshifts (GRB 090618 being one of 
the twelve in the list) to show that the time- 
resolved correlation (^E^^,^/^ — Ll^^^ within indi- 
vidual GRBs is consistent with time-integrated 
correlation (Ep^ak — Ltso) indicating that there 
must be a physical origin for the correlation and 
may not be an outcome of instrumental selection 
effects. 

3. DATA ANALYSIS 

We used Fermi /GBM and Swift /BAT data for our 
analysis. The GBM consists of twelve Sodium Io- 



dide (Nal: nx where x runs from to 11) and two 
Bismuth Germanate (BGO : by, where y runs from 
to 1) detectors. The Nal covers the energy range 
> 8 keV to 1 MeV while the BGO covers a wider 
range of ~ 200 keV to ~ 40 MeV. Data are saved 
in three types of packets (Meegan et. al. 2009): 
(1) TTE (Time Tag Events) in which time and 
energy information (in 128 energy channels) of in- 
dividual photons are stored, (2) CSPEC: binned 
data in 1.024 s bins and 128 energy channels and 
(3) CTIME: binned data in 0.064 s bins and 8 
energy channels (which is not suitable for spec- 
troscopy). The triggered detectors for the GRB 
090618 were n4, n7, bO, and bl. In the follow- 
ing analysis we have used n4, n7 and either bO or 
bl, one at a time. On the other hand, BAT has 
an array of CdZnTe detector modules located be- 
hind a Coded Aperture Mask. This covers an en- 
ergy range of 15 keV to 150 keV (Barthelmy et al. 
2005; Gehrels et al. 2004). The large area of BAT 
(5200 cm^) gives very good sensitivity for timing 
analysis, compared to that of Nal/BGO (126 cm^) 
detectors of GBM. 

3.1. Spectral Analysis 

Ghirlanda et al. (2010) have presented the time 
resolved spectral analysis results for GRB 090618 
in 14 time segments. We essentially followed the 
same methodology, but have also used the si- 
multaneous Swift /BAT data. GRB 090618 trig- 
gered the Swift /BAT instrument on 2009-06-18 
at 08:28:29.851 or 267006514.688 MET (s). The 
Fermi /GBM trigger time is To (GBM) = 2009-06- 
18 08:28:26.659 UT or 267006508.659 MET (s). 
Hence there is a delay time in the BAT data file 
which nmst be subtracted during the joint analy- 
sis. The calculated time delay is -3.192 s and when 
we convert it to MET it gives 267006511.496 MET 
(s) for Swift /BAT. In all the analysis this time is 
subtracted from the trigger time unless specified 
otherwise. 

For the GBM we chose the TTE (Time Tag Event) 
files of the triggered Nal detectors, namely n4 and 
n7. The BGO detector was excluded as it is not 
suitable for time resolved spectroscopy due to its 
lower sensitivity (Ghirlanda et al. 2010). The en- 
ergy range was chosen to be > 8 keV to ~ 900 
keV. We used RMFIT package provided by the 
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Fig. 1. — : Comparison of the parameters obtained from a time resolved spectral analysis of the prompt 
emission of GRB 090618 using a cut-off power law model. The left panel shows the power law indices 
obtained from a joint GBM and BAT analysis plotted against those obtained from a GBM analysis. The 
right panel shows a similar plot for the cutoff energy Eq. In both cases the straight line has a slope 1 denoting 
the equality of measurements (—a is plotted for convenience). 



User Contributions of Fermi Science Support Cen- 
ter {FSSC) to select the desired time cuts and for 
generating the background files. The background 
was taken before and after the burst, avoiding 
burst contamination (Ghirlanda et al. 2010) and 
the exposure time for background was optimized 
for different detectors. The energy channels were 
binned so as to get a minimum of 40 counts per 
channel. For the analysis of the BAT data we used 
standard tasks prescribed in The Swift BAT Soft- 
ware Guide, Version 6.3 of Heasoft 6.9 employing 
a time delay of -3.192 s. 

These spectra were then fitted with a cut-off power 
law of the form F = cons x E°'exp{—E/Eo) 
where F is the photon flux. The constant fac- 
tor, cons, was kept free while fitting essentially to 
take care of the different systematic errors in the 
area calibration of the two instruments. The de- 
rived parameter values (using GBM alone) were 
found to be consistent with the values presented 
in Ghirlanda et al. (2010). Figure 1 shows the 
comparison results of the spectral analysis. In Fig- 
ure 1 (left panel) the indices a obtained from the 
joint analysis are plotted against those obtained 
from an analysis of the GBM data for the four- 
teen different time segments. Similar data for the 
parameter Eq are shown in the right panel. A de- 
viation is noted for the parameter a in the to 3 
s region. This may be attributed to lower count 



rate of BAT in the precursor region. 

We note a systematic increment in the value of 
a and Eq after the inclusion of BAT data. Fur- 
ther, the constant factor is always lower by 10-20% 
for the BAT. Sakamoto et al. (2011) have made a 
joint spectral analysis of 14 GRBs using data from 
Swift/BAT, Konus/Wmd and Suzaku/WAM and 
they have noted that for the BAT data, the con- 
stant factor is systematically lower by 10 - 20%, 
the photon index is steeper by 0.1 0.2 and the 
Epeak becomes higher by 10 - 20%. The results 
presented here are in general agreement with these 
conclusions. 

In Figure 2 we have plotted the ratio of fractional 

errors as a function of the measured parameters. 
We find that there is a marginal improvement for 
the parameter Eq (Figure 2, right panel). The pa- 
rameter a (Figure 2, left panel), however is mea- 
sured always with a greater accuracy from the 
joint analysis. The large effective area of BAT im- 
proves the accuracy of the spectral index whereas 
its low high energy threshold (~150 keV) makes 
relatively low impact on the measurement of the 
high energy cut-off value. The joint analysis al- 
ways gives confidence on the measured parameters 
and ensures that the effects of instrumental arte- 
facts are minimized. Hence we use, in the later 
analysis, a joint fit for measuring the global spec- 
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Fig. 2. — : The ratio of fractional errors for GBM-BAT joint analysis to the GBM analysis is shown as a 
function of the measured values of a (left panel) and Eq (right panel) obtained from a fitting for GBM data 
alone. The dotted line in both cases shows the average of all errors measured. This is found to be marginally 
less than the line of equal error (dot dashed) for the Eq and well below that line for a measurement (—a is 
plotted for convenience). 



tral parameters and use only the BAT data for a 
joint timing and spectral analysis. 

3.2. Timing Analysis 

We identified four pulses (other than the precur- 
sor) in the ^60 s to ^130 s period after the trigger. 
Both the Swift /BKT and Fermi /GBM (n4) were 
used for the timing analysis. We subtracted back- 
ground and shifted the time origin to the trigger 
time (TRIGTIME is the keyword in the header) 
of the respective satellite. As we noted previously, 
there is a time lag of -3.192 s for the BAT trigger. 
Light curves were extracted for different energy 
bands for both the data sets. The energy bands 
for the Swift /BAT were: 15-25 keV, 25-50 keV, 50- 
100 keV, and 100-200 keV. For the Fermi /GBM 
(n4) we chose 8-15 keV, 15-25 keV, 25-50 keV, 50- 
100 keV, 100-200 keV, and 200-500 keV. Greater 
than 500 keV region of Nal and the BGO detector 
of GBM were excluded because the third and the 
fourth pulses were not detected in these bands. 

All the light ctirves were fitted with the four- 
parameter model of Norris et al. (2005) 

I{t) = AnXexp{-n/{t - t,) - (i - ts)/T2} (1) 

for t > ts, where = (Ti/r2)^ and A = exp{2(x). 
An is the pulse amplitude, ts is the pulse start time 



and Ti, T2 are time constants characterizing the 
rise and decay part of the pulse. From these pa- 
rameters, various derived quantities like the peak 
position, the pulse width (measured as the sepa- 
ration between the two 1/e intensity points), and 
the pulse asymmetry can be calculated (see Norris 
et al. 2005 for details). Errors are calculated for 
the parameters as nominal 90% confidence level 
errors (Ax^=2.7) and these errors are propagated 
for the derived quantities, as described in Norris 
et al. (2005). 

Figure 3 shows the pulse width (w) as a function 
of energy (E) for all the pulses. Apart from the 
fact that the width broadens in the lower energy 
bands, we note a tentative evidence for a decre- 
ment (or constancy within error bar) of the width 
from 50-100 keV through 8-15 keV energy bands 
for the third and the fourth pulses (lower panels 
of Figure 3). We fitted these data for the individ- 
ual pulses with (i) constant {w = c) and (ii) linear 
{w = mE -\- c) models. For the first two pulses 
there was a significant improvement in the for a 
linear fit (the varied from 16.2(9) to 5.1(8) and 
140.9(9) to 7.9(8) respectively for the two pulses; 
the numbers in the brackets are the degrees of free- 
dom). The derived slopes are (—7.4 ± 2.6) x 10~^ 
s keV~^ and (-5.5 ±0.4) x 10~^ s keV"\ respec- 
tively for the two pulses. The third and the fourth 
pulses, however, showed negligible improvement in 



5 



8 

7 

6 

20 

15 

10 
L 

15 
10 











H 1 1 1 1 1 1 1 1 1 1 1 ^ 

+ ~ * 

1 1 1 1 1 1 1 1 1 1 1 1 


H 1 ^ 




'-^^ : t : (, - 



100 



200 



300 




50 100 150 

Energy (keV) 

Fig. 3. — : Pulse width (w) is plotted against en- 
ergy (E), for pulse 1 to 4 respectively in the top to 
bottom panels. The fourth pulse is insignificant in 
the 200-500 keV energy range. Hence the fourth 
panel has a different scale on the x-axis. Values de- 
rived from the Swift /BAT and the Fermi/GBM 
data arc shown by circles and filled circles, respec- 
tively. While the width decreases with increasing 
energy for the first two pulses (top two panels), 
there is a tentative evidence of the reverse effect 
in the lower energy (~ 10-70 keV band) for the 
other two pulses (see section 3.2). 



lection of integration time for pulse spectral anal- 
ysis of the individual pulses. The pulse spectral 
analysis was done with a Band function fit and in- 
cluding the BGO detector of GBM. We included 
the precursor pulse in this analysis. We divided 
the GRB event into four time segments namely. 
Part 1 (To to Tq + 50), Part 2 (To + 50 to To + 77), 
Part 3 (To + 77 to Tq + 100) and Part 4 (Tq -t- 100 to 
To -|- 180). The first and second pulses are consid- 
ered together as there is too much overlap between 
them. 

All the spectra were fitted in the energy range ^ 
8 keV - ~ 1 MeV with the four parameter model 
of Band (Band et al. 1993): 



m = { 



if £ < [(a - /3)/(2 + a)]Ep,ak 



Ab [ exp [0 - a] 



(o!-^)-Epeafc 

100(2-|-a) 



otherwise 



(a-/3) 



(2) 



(from 8.7(9) to 6.0(8) and 2.1(8) to 2.1(7), re- 
spectively for the two pulses) for the linear fit. For 
these two pulses, when the data below 70 kcV is 
considered, there is a marginal improvement in the 

for a linear fit with a positive slope (from 8.7(9) 
to 3.8(5) and 2.1(8) to 1.5(5), respectively for the 
two pulses). The derived slopes are (15±14) x 10~^ 
s keV-i and (10 ±21) x IQ-^ s keV^^ respectively 
Though the slopes for these two pulses are posi- 
tive, the error bars are quite large and hence the 
evidence of this reverse variation is tentative. In 
Figure 3 (lower two panels) the linear fit for pulse 
3 and 4 are shown by solid lines (in the low en- 
ergy) and dot-dashed lines (in full energy band). 
We have ignored the fourth pulse in the 200-500 
keV band of Fermi for the fitting as it was not 
prominent in this energy band. 

3.3. Pulse- Wise Joint Spectral Analysis 



The lower and the higher energy parts are char- 
acterized by the parameters a and /3 respectively 
and they meet smoothly at the spectral break en- 
ergy [(a - /3)/(2 + a)]Epeak- The peak energy, 
Epeak or equivalently the characteristic energy, de- 
fined as Eq = Epeak/ (2 + a), may be considered 
as the third parameter of the model. The best 
fit model parameters along with the nominal 90% 
confidence level error (Ax^=2.7) for the joint fit 
of Fermi/GBM (n4, n7, bl) and Swift /BAT are 
reported in Table 1 (bO in place of bl gives sim- 
ilar results). The high energy photon index /3 is 
constrained primarily by the GBM data and hence 
we have kept its value fixed at -2.5 to determine 
precise values of Ep^ak and a. Only in one case 
namely for the full range (Tq to Tq + 180 sec) of 
the GRB, /3 was made to vary and a value of - 
2.96 was obtained. Figure 4 shows the joint Band 
spectral fitting of this range along with the de- 
viation of the data. 



The peak position and the width gave a clue for 
the time cut for each pulse and facilitated the se- 



An examination of the residuals in Figure 4 reveals 
that there are some systematic differences in the 
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Table 1:: spectral parameters of GRB 090618 for Band model, from a joint analysis of Fermi/GBM (n4, n7, 
bl) and Swift /BAT data 



Part 



/3 



^peafc(keV) x"(dof) Eq (keV) 



Full 


-1.35±0.02 


-2.5 (fixed) 


154±7 


428(299) 


236±18 


Part 1 (Precursor) 


-1.25±0.05 


-2.5 (fixed) 


166tJ| 


369(299) 


221^39 


Part 2 (Pulse 1 & 2) 


-1.11±0.02 


-2. 5 (fixed) 


212±9 


460(299) 


238±15 


Part 3 (Pulse 3) 


-1.15±0.03 


-2.5(fixed) 


114±4 


334(259) 


134±9 


Part 4 (Pulse 4) 


1.00_o.o8 


-2.5 (fixed) 


33±3 


431(299) 


04+38 


Full (free p) 


-1.36±0.02 


-2.96±0.48 


160t^ 


418(288) 


250^22 



response in the 20-40 kcV region between BAT and 
GBM data and the resultant 




20 50 100 200 500 

Energy (keV) 



Fig. 4. — : Fitting of Band function to the full 
range (0 to 180 s) of GRB 090618 for the energy 
range ~ 8 keV to ~ 1 MeV. The range of the BGO 
detector starts above 200 keV, while BAT covers 
up to ^ 150 keV. The other two detectors namely 
n4 and n7 have the spectrum above 8 keV to ~ 
800 kcV. The extreme channels and specifically 
the 25th channel of n4 are ignored due to high 
background rates. The unfolded spectrum of BAT 
shows a lower normalization (see text). 



values generally are quite high. The resultant 
spectral parameters, however, agree with the re- 
sults reported for this GRB (Ghirlanda et al., 
2010; Nava ct al. 2011; Rao et al. 2011). The er- 
rors are significantly better compared to the joint 
BAT and RT-2 analysis (Rao et al. 2011), but 
comparable to that obtained from a GBM analy- 
sis alone (Nava et al. 2011). Since a joint analysis 
using data from different satellites pins down the 
systematic errors in each instruments, we believe 



that the joint analysis results, reported here, are 
devoid of instrumental artefacts, particularly for 
the parameters a and Ep^ak- 

4. JOINT SPECTRAL AND TIMING 
STUDY: AN ALTERNATE APPROACH 

4.1. Two Parameter XSPEC Table Model 

4.1.1. The time evolution of Ep^ak 

The light curve of each pulse of a GRB can be 

described in terms of the four-parameter model 
of Norris (equation 1) and the spectrum for each 
pulse can be fitted with the Band model (equation 
2). Liang and Kargatis (1996) showed for clean, 
separable pulses that the peak energy of the vF^ 
spectrum, Ep^ak, decays with the photon fluence 
as, 

Epeakit') = EpeakM e-*^''^/'''° (3) 

where Ep^akfl is the peak energy at zero fluence, 0o 
is the decay constant and t' = t — tg- The quantity 
4){t') is the fluence up to time t from the start of 
the pulse {tg). For convenience of description we 
shift our origin to tg and hence t' — t for each 
pulse analysis. Equation 3 quantifies the overall 
softening of the burst with time. 

Kocevski and Liang (2003) have integrated the 
light curves of a sample of GRBs up to various 
times t and found the fiuence, 0(i), at those times. 
The peak energy at the same times, Ep^akit), 
were found by fitting the spectra with a Band 
model. They determined the Epeak,o and (po from 
a semi-log plot of Epeak{t) versus (j)(t) and fitting 
a straight line. This procedure was adopted for 
clean, bright, and separable FRED pulses. 
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Here we explore a method to determine the timing 
and spectral parameters by assuming that Epeak{t) 
and 4'{t) are related according to equation 3 in each 
pulse. Essentially we assume that for each pulse, 
the variation of flux with time is given by the ex- 
ponential model (equation 1), and at each time the 
spectra is described by a Band function (equation 
2), and, for a given pulse, only the peak energy 
varies with time governed by equation 3. We de- 
velop a method to derive the model parameters 
from a time integrated spectrum and energy inte- 
grated light curve. 

4.. 1.2. Description of the method 

The flux of a GRB pulse is a function of time and 
energy. If there exist empirical formulae for the 

temporal and spectral beahvior then the flux can 
be described in terms of all the model parame- 
ters. For example, the variation of flux with time 
is given by Norris' exponential model, with four 
parameters (A„, the normalization, tg^ ti, and T2). 
The spectral behaviour is parametrized by a Band 
function, with four parameters {A}j, the normaliza- 
tion, a, /?, and Epeak)- For any timing model, the 
time parameters change with energy e.g., the de- 
rived parameter, width, changes in various energy 
bands (see section 3.2). Similarly, the spectral pa- 
rameters of any spectral model change with time 
e.g., a and Eq of cut-off power law change with 
time in the same pulse (see section 3.1). 

The Norris exponential model (equation 1) de- 
scribes the evolution of a pulse with time. This 
description can be used for counts in any energy 
band and the variation of the model parameters 
with energy can be studied. Instead, we assume 
that the model parameters ti and T2 derived for a 
pulse integrated over all energy bands are global 
parameters describing the evolution of total flux 
with time. Since the variation of spectrum is best 
described by the peak energy evolution with flu- 
ence (equation 3), we further assume that a and 
/3 do not change with time. The normalization 
factors {An and Ab) are related by the fact that 
the fluence calculated either by time integration 
of the pulse light curve or by the integration over 
the total energy should be the same. With these 
assumptions, the parameters for a given pulse are 
determined as follows. 



For a given pulse, we use the energy integrated 
light curve from BAT and use equation 1 for the 
fitting. This best fit light curve is used to deter- 
mine the global parameters namely. An, ti, and 
T2. Since we can arbitrarily choose the pulse start- 
ing time by shifting the time origin, the parame- 
ter ts is not very relevant (also see section 4.1.1). 
This function is then integrated up to a time t to 
determine (p{t). From the individual spectral fit- 
ting for each pulse using the Band model, a and 
(3 are determined (see section 3.3). Now, if the 
constant parameters namely, Epgakfl and (po are 
determined, we have a complete spectral and tim- 
ing description of a GRB pulse. For this purpose, 
we use a two-dimensional grid of guess values for 
these two parameters and then determine the best 
fit values by a spectral fitting method. 

For a particular combination of i?peafe,o and (f>Q 
values, Epeak{i) at any time t can be determined 
from equation 3. On the other hand, equation 2 
can be rewritten as I{E) = Ai, f{E) where Ai, is 
the normalization and f(E) is the rest of the func- 
tional form valid in the two energy bands. Given 
the global parameters {a and /3) and the time de- 
pendent parameter (Epeakit)), f(E) can be deter- 
mined. The normalization Ai,{t) can be found by 
dividing the fluence (</>(t)) of the light curve by 
the f(E) integrated over the whole energy range. 
Hence photon spectrum is generated for the time 
instance t. Choosing a time and energy grid, 
counts can be plotted as a function of time and 
energy for a particular combination of Ep^ak.o and 
00 values, and the resulting surface can be referred 
to as a three dimensional pulse model. Hence for a 
two-dimensional grid of guess values of Epeak.o and 
(j)Q we can generate three dimensional pulse mod- 
els for each of the grid points. These model pulses 
can be used in two ways: (1) integration over time 
gives the model spectrum of that pulse, (2) in- 
tegration over various energy bands gives model 
light curves of that pulse in those energy bands. 

In our analysis, we first generate a two param- 
eter XSPEC table model using the guess values 
of Epeakfl and <j)o and then use this table model 
to determine the best fit values of these two pa- 
rameters by fitting the time integrated spectrum 
(see section 4.1.3). With this method, we get a 
complete description of the pulse with model pa- 
rameters determined from the data. 
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Table 2:: Best fit values of Ep^akfi and for the pulses as determined by the minimization of the 
spectral data with the additive table model. The first two pulses are combined together. The other pulse 
characteristics (n, T2, and ts) are also reported. 



Pulse 






Norm 


X' (dof) 


T1(S) 


T2(S) 


i.(s) 


1 


35911^ 




0.74 ±0.06 


40.12 (75) 


795.4+^-^ 




4o.it;^-i 


2 










858.2+?;^ 


n t:o+U.UU4 
^•^"-0.005 


^^•^-0.2 


3 


324;«^ 


18.8+^2 


0.54 ±0.04 


31.41 (75) 


353.4+*g 


2 47+'^-i^ 

^•^'-0.09 


50.3+i-^ 


4 


307t- 




0.19 ±0.04 


56.42 (75) 


532.0i^^^ 







'4 




Fig. 5. — : The three dimensional model of the 

whole GRB using a model with the best fit values 
of Epeakfl and (j)o for each pulse. Counts in ar- 
bitrary unit is plotted in the vertical axis against 
time and energy. The pulses are shifted to original 
pulse start time (see Table 2). 



Using this model we can reconstruct the three 
dimensional pulse profiles and can generate light 
curves in any energy band. Various timing param- 
eters (width and lag) are derived in these energy 
bands (see section 4.2) and these model predicted 
values can be compared with the observed values 
to get a validation of the whole procedure. 

4.1.3. Specification of the table model and the 
best fit values 

For all the pulses we use a time grid resolution of 
0.5 s and energy grid resolution of 2.0 keV. Each 
pulse has a separate XSPEC table model, though 
the first and second pulses are combined together 



as the data overlap. A typical table model con- 
tains a spectrum in the 2-200 keV energy range for 
200 values of Epeakfi from 100.0 keV to 1100.0 keV 
and 50 values of ^0 from 2.0 photon cm~^ to 77.0 
photon cm~^, energy binning information, param- 
eter values and the required header keywords stan- 
dardized as per the requirement by XSPEC de- 
scribed in OGIP Memo OGIP/92-009. 

The additive table model is used to determine the 

best fit values of Ep^akfi and (j)^ by the standard 
minimization procedure. The best fit values 
with nominal 90% confidence errors are reported 
in Table 2. The parameters, Epeak.o and 4>o of the 
first two pulses are assumed to be the same. It 
should be noted here that we have used the nor- 
malization A global parameter to calculate 
4>{t). The normalization of the subsequent pulse 
profile is derived from the integration of f(E) over 
the total energy band (see section 4.1.2). Effec- 
tively this leads to a spectrum with normalization 
as a free parameter. This normalization is deter- 
mined during the spectral fitting using XSPEC, 
and is also reported in Table 2. 

4.1.4. Three dimensional pulse model for the best 
fit values 

Three dimensional model of each pulse is gener- 
ated for the best fit parameter values (as described 
in sections 4.1.2 and 4.1.3). The three dimensional 
model of the whole GRB obtained from the sum of 
these models of the pulses with each pulse shifted 
to its proper starting time is shown in Figure 5. 
This figure is used only for display purposes. In 
actual analysis, each pulse is normalized properly, 
before calculating various parameters as the over- 
all normalization may vary for each pulse. This 
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Fig. 6. — : Observed light curve in the full energy range 15-200 keV fitted with the Norris' exponential model 
is shown in the left panel. The timing parameters are determined and used in the pulse generation to predict 
light curves in all energy bands. One such predicted light curve in the 25-50 keV region is shown along with 
the observed light curve in the right panel. The count/s in the y-axis is BAT mask weighted count s— 1. 



normalization is achieved for spectrum and light 
curve separately. 

4.2. Timing analysis using the model 

4-2.1. Reproducing the light curves 

Once we get the best fit values of Epeak.o and cpQ for 
each pulse we can generate the best fit model light 
curve by integrating the three dimensional pulse 
model for the best fit value over the desired energy 
range. We generate three dimensional pulse mod- 
els, shift them over appropriate starting time, co- 
add them with proper normalization factors and 
integrate over energy ranges to generate the com- 
bined light curves in the energy bands: 15-25 keV, 
25-50 keV, 50-100 keV and 100-200 keV of BAT. 
In Figure 6 (right panel), we show the observed 
light curve along with the model predicted light 
curve for the 25-50 keV energy band. In the left 
panel Norris' exponential model fit for the total 
energy range is shown. Parameters from this fit 
(Figure 6, left panel) were used for characterizing 
the timing parameters in our model. 

The GRB light curve shows rapid time variable 
components over and above the smooth pulse 
structure. In the pulse description we consider 
only this smooth variation. Any small rapid time 



variability in the data gives large value of when 
we use the smooth pulse model for fitting. Hence, 
the normalization factor could not be determined 
by minimization. Instead we physically ex- 
amined the heights and determined the normal- 
ization. The co-added pulses reproduce the light 
curve quite well. Hence, we can conclude that the 
total light curve and the integrated spectrum, with 
the assumption of £^peafe(i) varying with <p{t) in a 
predetermined way (equation (3)), correctly repro- 
duces the energy resolved light curves. This can be 
taken as a confirmation of the assumption of the 
Epeak variation. It should, however, be noted that 
a single joint spectral fit is done for pulses 1 and 2 
and hence the derived values of i^peafe.o and 0o for 
the first two pulses are only the average values. 

4.2.2. Derived timing parameters 

Pulse width of individual pulses in the various en- 
ergy bands can be derived by fitting the light curve 
using equation (1) and using the best fit values of 
Ti and T2 (see section 3.2). For overlapping pulses 
this method is erroneous e.g., rising part of one 
pulse can fall on top of the falling part of the pre- 
vious pulse affecting the width. Also, pulse delay 
characteristics are affected by these overlaps. In 
such cases the method described here is very use- 
ful. It takes as input some global parameters {An , 
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Tabic 3:: Variation of the model predicted width of each pulse with energy. Numbers within brackets are 

the measured pulse widths. 



Pulse 


15-25 keV 


25-50 keV 


50-100 keV 


100-200 keV 


1 


(6.74t°;}|) 


(6.67t°;f9) 


(6.i8l};i) 


5.08t;;:i^ 

(5.7910:30) 


2 


6.56115:- 

(7.52l°;i^) 


6.35i^;^^ 


5 89+'^-^''' 

'->-o»_0.34 

(7.23l°;il) 


r, 00+0.20 
^••^"^-0.34 

(6.821°:°^) 


3 




is.io^U.i 
(i8.09i°:^0) 


18.68+-- 
(17.22tl:^^) 


(15.87ll:i) 


4 


(13.32t°jJ) 


14 66+'^-^^ 

i'i.OO_Q gg 

(13.76tlJ0) 


10.iO_Q gg 

(14.0011:23) 


15.93+i-^? 



Ti, T2, a and /?) which are average quantities char- 
acterizing the pulse and some constant parameters 
(Epeak.o, '/'o), dctcrmincs the variable parameter 
{Epeak{t)) over the time grids and generates the 
three dimensional model of the pulse. Note that 
the timing parameters ti and T2 are determined 
from the energy integrated light curve and are 
fixed in our model. The timing properties of the 
pulse (e.g., width in various energy bands, spectral 
delay) are now determined crucially by the evolu- 
tion of the peak energy. Integration of the three 
dimensional model in various energy bands gener- 
ates the light curves in those energy bands. Hence 
physically we can determine the full width at the 
1/e intensity and the peak positions. Wc define 
the spectral delay as the shift in the peak posi- 
tion at each energy band compared to the lowest 
energy band (15-25 keV). 

The derived pulse width and delay are reported in 

Table 3 and Table 4, respectively. In Table 3 the 
numbers in the brackets are the width determined 
from Norris' exponential fitting for individual light 
curves. Figure 7 shows a comparison between the 
measured width following the new method with 
that determined by fitting the light curve. We note 
a systematic decrement of the width measurement 
by our method for the first two pulses. This is 
simply due to the fact that individual pulses rather 
than co-added pulses are taken into consideration. 
Hence, measured widths are more accurate and 
devoid of contamination effects from other pulses. 

One of the pulse properties of GRB is broaden- 
ing at lower energies (Norris et al. 1996, Hakkila 
2011). Figure 8 (left panel) shows the pulse width 



25 



20 - 




5 10 15 20 25 

Width (New method) 

Fig. 7. — : Comparison of the width of individ- 
ual pulses in various energy bands determined by 

directly fitting the light curves (Norris model fit) 
and the three dimensional model prediction (new 
method). 



variation with energy. While the first two pulses 
follow the normal trend, we notice that the width 
for the third and fourth pulses broadens with 
higher energy. This cannot arise from the contam- 
ination of the previous pulse as could have been 
the argument for fitting the entire light curve. In 
our case, each pulse is generated and analyzed sep- 
arately for the width determination. The energy- 
width plot of individual pulses is fitted with a 
linear model. The slopes are negative for the 
first two pulses ((-9.4 ± 1.7) x lO'^ s keV'^ and 
(-11.8 ± 2.4) X 10^3 s keV-i respectively; = 
0.13 (2) and 0.11 (2) respectively) and positive for 
the other two ((18.8 ± 11.1) x 10~3 s kcV'^ and 
(14.2 ± 10.6) x 10-3 s keV-i respectively; = 
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Tabic 4:: Variation of the model predicted delay with energy for each pulse. Delay is calculated with respect 
to the 15-25 keV band (mean energy: 20.82 keV). 



Pulse 


Energy Channel 


Mean Energy (keV) 


Delay (s) 


1 


15-25 vs. 25-50 keV 


35.45 


0il35 n r\oo 
— U.UZZ 




15-25 vs. 50-100 keV 

15-25 vs. 100-200 keV 


68.07 

123.73 


n 07C+0.023 

_Q 7OA + O.O43 


2 


15-25 vs. 25-50 keV 


35.45 


-U.iyti_Q 




15-25 vs. 50-100 keV 


68.07 






15-25 vs. 100-200 kcV 


123.73 




3 


15-25 vs. 25-50 kcV 


35.45 






15-25 vs. 50-100 keV 


68.07 


-0 04n+" "i2 

U.U4U_Q 030 




15-25 vs. 100-200 keV 


123.73 


-n nSr,+0.005 
U.UOO_o.005 


4 


15-25 vs. 25-50 keV 


35.45 


-0 n20+"-o^u 

u.uzu_o 010 




15-25 vs. 50-100 keV 
15-25 vs. 100-200 keV 


68.07 
123.73 


-0 065+°-0^5 

U.UDO_o 

_n nqc;+0 005 
u.uyo_o.oo5 



0.0046 (2) and 0.0021 (2) respectively). Figure 
8 (right panel) shows a comparison of the slopes 
meaurcd by the linear fit of data determined by 
the two methods (compare slopes determined in 
section 3.2). We note again that the slopes of the 
first two pulses are convincingly negative (low er- 
ror bar), but slope of the other two pulses have 
large error bars. Hence, the evidence of reverse 
pulse broadening is tentative, though individual 
pulse analysis in the new method shows a bet- 
ter evidence for the reverse pulse broadening effect 
(see Figure 8 - left panel). 

This deviation from the canonical pulse width 
broadening with energy decrement may have a 
physical significance. In our analysis physical pro- 
cess of pulse emission is not considered. But 
phenomenologically, this occurs depending on the 
values of ri and ti. A detailed comparison be- 
tween the first and the third pulse reveals the fact 
that they have different pulse width variation al- 
beit having nearly similar parameters except for t\ 
and Ti. These values are as follows (the numbers 
within parentheses are those of the third pulse): 
a= -1.11 (-1.15), /3 = -2.5 (-2.5), i;peafe.o=359 
(324), 00=12.2 (18.8), ri=795.4 (353.4), ti =0.54 
(2.47). It is apparent that the main contributing 
factors for the different width-energy relation are 
the global values of t\ and ti . Hence we examined 
for a generic pulse keeping all other parameters 
fixed at the values of the first pulse and varying n 
for different values of ti. 



Figure 9 shows the regions in which pulse broaden- 
ing or its reverse phenomenon can occur. This fig- 
ure is generated by fixing all the parameters other 
than n and ti to those for the first pulse. For four 
values of T2 wo vary t\ and note the ratio of the 
width in 100-200 kcV region (104) to that in the 
15-25 keV {w\) region. The solid line parallel to 
Ti with W4/wi=l divides the plot into two regions: 
(i) normal width broadening region for which Wi 
< wi and (ii) anomalous width broadening region 
for which W4 > Wi. Global values of ri and T2 
falling in the W4 < Wi region will lead to normal 
pulse broadening (i.e., pulse broadens with lower 
energy) and vice versa. 

We identify from the figure that the ti and T2 val- 
ues (795.4 and 0.54 respectively) of the first pulse 
lie well within the region of normal width broad- 
ening. We also note that the values of ti and T2 
for the third pulse (353.4 and 2.47 respectively) lie 
within the region of anomalous width broadening 
despite the fact that the parameter values used, 
other than n and T2, are those of the first pulse. 
This shows how insensitive is the pulse broaden- 
ing effect with all other parameters and crucially 
determined by the global values of ri and T2. 

Another pulse property is the hard-to-soft spec- 
tral evolution. In Figure 10 we have plotted the 
model predicted delay as a function of energy (left 
panel). It can be seen that the spectral time delay 
always shows a decrement with energy. The hard 
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Fig. 8. — : (Left panel) Pulse width (w) variation with energy (E) as predicted by our model. Pulse 1 to 4 
are shown from the top to bottom panels. Width broadens with lower energy for the first two pulses (normal 
width broadening) while the reverse effect is apparent for the third and the fourth pulses (anomalous width 
broadening). (Right panel) Comparison of slopes of E-w curves determined by Norris model and the new 
method. 



X ray always precedes the softer one. The val- 
ues for the first three peaks can be compared with 
the observed delay given in Rao et al. (2011) and 
this is shown in the right panel of Figure 10. The 
time intervals taken by Rao et al. (2011) are Tq to 
To+50 s, To+50 s to Tq+T? s, To+77 s to To+lOO 
s, To+lOO s to To+180 s. Among these the first 
time bin contains the precursor burst which is not 
taken in our calculation. The first and the second 
pulses are combined in the second time bin. The 
third pulse is covered in the third time bin. Data 
of the fourth pulse cannot be compared with the 
time bin as they differ by ~ 55 s (80 s in time 
bill as compared to ~ 25 s in pulse analysis). The 
deviation of the model predicted delay from the 
observed delay can be explained by the fact that 
the delay in the second time bin (To+50 to To+77) 
is the combined effect of the first two pulses. The 
data for the third time bin agrees within error with 
the model predicted delay. Also, the observed de- 
lays are calculated using peak position measure- 
ment and cross correlation of the pulses whereas 
model predicted delays are calculated from peak 
position measurement alone. 



4.3. Correlations and the Bigger Picture 




2000 



Fig. 9.—: Ratio of the width in 100-200 keV {W4) 
to the width in 15-25 keV (wi), as predicted by 
the pulse model, is plotted against ti for different 
values of T2. In the figure maximum as well as 
the minimum of the curves move to lower ti val- 
ues for higher T2. The values of T2 are 0.54, 1.20, 
1.90, and 2.47, respectively from left to right. The 
Ti values determined for the second and the third 
pulses of the GRB (795.4 s and 353.4 s, respec- 
tively) are shown by dash-dotted lines. The solid 
line divides the region into two parts: the Nor- 
mal width broadening region and the Anomalous 
width broadening region (see section 4.2.2). 



One of the important properties of the prompt 
emission is that the various measured quantities 



correlate with luminosity. Particularly, the peak 
energy of the vFi, spectrum, Epeak correlates with 
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Fig. 10. — : Model predicted delay of light curve with mean energy E, with respect to the 15-25 keV light 
curve (left panel). For the four pulses the data are shown in square, circle, triangle, and star respectively. 
Right panel shows the comparison of the measured delay (Rao et al. 2011) with the model predicted delay. 
The lower six points are model prediction for first and second pulses compared with the 50-77 sec interval 
(see text). The end points along the x-axis for this data group are the model predicted delays with the left 
one for the second pulse. The measured delay along with error is shown in the y axis. The straight line is 
the line of equality. 



isotropic energy Eis,, (Amati et al. 2002) and also 
with the peak isotropic luminosity Lisa (Yonetoku 
et al. 2004). Ghirlanda ct al. (2004) showed 
that the correlation is tighter with collimation- 
corrected energy, E^^. Schaefer (2007) used a sot 
of 69 GRB sample (up till 7 June 2006) with 
measured redshift to calculate various correla- 
tions viz. lag- luminosity {riag — L), variability- 
luminosity {V - L), Epeak - L, Epeak - E^, min- 
imum rise time-luminosity {trt — L), number of 
peaks-luminosity (Npeak — L) to get distance mod- 
uli for each. Weighted average of them gave the 
distance modulus {fj) which is plotted against the 
measured redshift (z) extending the Hubble Dia- 
gram (HD) to z >6. Therefore, these correlations 
can be used to standardize the GRB energy budget 
making GRBs a class of cosmic distance indica- 
tors. Also, these correlations might point to a un- 
derlying fundamental physical process (Yamazaki 
et al. 2004; Rees & Meszaros 2005; Xu et al. 2005; 
Firmani et al. 2005; 2006; 2007; Ghirlanda et 
al. 2006a;b; Wang et al. 2006; Thompson 2006; 
Thompson et al. 2007; Liang et al. 2005; 2008 
and the references therein; Li et al. 2008; Qi et al. 
2008). 

In order to establish the reality of these correla- 



tions one must show that all of these correlations 
hold just as the same way within a single GRB 
as for a sample of many GRBs. Ghirlanda et al. 
(2010) showed for a set of GRBs that the time- 
resolved correlations of a single GRB are similar 
to the time integrated correlations among differ- 
ent GRBs. In the present work we have used the 
underlying pulse structure to investigate the cor- 
relations: Epgakfi ~ Liso, Epgak,0~ Ej^igo, 4>0~Liso 

and (/>o — E^^iso within the pulses of a single GRB. 

The correlations are shown in Figure 11. The 
mean and the Su scatter of Epeak — Lisa (top 
left panel) and Epeak — E^jso (top right panel: 
Ghirlanda et al. 2010) are shown for reference. 
As expected, the Epeak,o values always lie above 
the mean of Epeak values. It seems that the corre- 
lations of Epeakfi with Ljgo and E^^iso are rather 
tighter in our case, though with only three points 
we cannot draw any definite conclusions. While 
the parameter Epeak,o is a constant, Epeak is an av- 
erage quantity. Hence, we expect that Epeak,o Pro- 
vides a better pulse description than Epeak- This 
can be tested in future for correlations of Epeakfi 
rather than Epeak with other quantities in differ- 
ent GRBs. But if the same correlations hold for 
different GRBs then we might conclude (a) the 
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Fig. 11. — : Correlations between {B)Epeakfi — Lisa (left panel, top); the dot-dashed line is the mean iJpeafe — 
Liso curve while the dashed lines represent 3cr scatter (from Ghirlanda et al. 2010), — Li^o (bottom); the 
dot-daslicd line is the mean i/io — Liso line (from Kocevski et al. 2003) and (b) Ep^akfi — E^,iso (right panel, 
top); the dot-dashed line is the mean Epeak — E^^iso curve while the dashed lines represent 3a scatter (from 



Ghirlanda et al. 2010), 



E-y^iso (bottom). 



correlations are rather physical and not an out- 
come of the bias in sample selection, (b) in a single 
multi-peaked GRB each pulse is generated inde- 
pendently. We also note that the calculated Liso 
and E^_iso are somewhat lower than those calcu- 
lated by Ghirlanda et al. (2010). This must have 
come from the instrumental effect as noticed ear- 
lier. The normalization of BAT is always less than 
that of Fermi. This leads to a lower flux calculated 
from the BAT spectrum. 

5. DISCUSSION AND CONCLUSIONS 

We have attempted to explain the timing and 
spectral variations within a pulse of GRB by a 
simple empirical model and have given a prescrip- 
tion to measure the parameters of this model. 
We use simultaneous data from Swift /BAT and 
Fermi/GJiM to determine the average pulse spec- 
tral characteristics using the Band model. Then 
we use the integrated light curve from Swift /BAT 
to get a simple description of the time evolution 
using Norris' exponential model. Using these pa- 
rameters, we develop a method to fit the spec- 
trum using XSPEC table model by assuming that 
Epeak{t) decreases exponentially with fluence in 
the pulse. We determine the parameters of the 
model and then find that these parameters can 



correctly describe the energy dependence of the 
pulse characteristics like width and spectral lags. 
We also find, in the data, a tentative evidence of 
anomalous width broadening with energy for some 
pulses. The joint spectral and timing analysis 
developed here confirms this phenomena for the 
same pulses at least in the lower energy bands. 

The method developed here is applicable for both 
single pulse and structured GRBs. Most of the 
long GRBs are structured which makes them dif- 
ficult for pulse analysis. As far as the correlations 
are concerned, pulse analysis is more indicative 
than intensity guided time-resolved spectroscopy. 
Hence to draw conclusion regarding these corre- 
lations within a sample of GRBs and within the 
pulses of a single GRB we must include these 
structured GRBs with multiple pulses. Kocevski 
et al. (2003) have prescribed a method for short, 
clean pulses. Our method extends the domain to 
all kinds of GRBs. Also the present method is 
much more physical than the time-resolved one in 
the sense that the individual pulses are treated 
separately. Hence there won't be contamination 
effects. The Ljgo-Epeafc.o looks tighter because 
correlations within pulses may be more physical 
than the time-resolved correlation or the time in- 
tegrated correlations. In future the method devel- 
oped here can be used for a large set of GRBs (onto 
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their individual pulses) to investigate the improve- 
ment in the correlations. 

It is very interesting to note that the two distinct 
pulses (pulse 3 and pulse 4) in GRB 090618 have 
Ejso values within a factor of 2, but widely dif- 
ferent Epea/c values (114 keV and 33 keV, respec- 
tively). They will show a huge scatter if used in- 
dividually for a global correlation. The ^peak,o 
values, on the other hand, are closer to each other 
and lead to a lower scatter in the correlation (see 
Figure 11. top right panel). Hence we suggest that 
the empirical relations are intrinsic to the emission 
mechanism of GRBs and Epeafc,o values, instead of 
Epeafe , should be the correct parameter to be used 
for correlation studies. 

The method described in this paper is quite 
generic in the sense that any empirical descrip- 
tion of the spectrum and its evolution with time 
can be tested against the data and the derived pa- 
rameters can be used for global correlation studies. 
For example the Band spectrum is derived from 
an extensive study of data from CG'i?0/BATSE, 
which is primarily sensitive above 50 keV (Band et 
al. 1993). Extension of the data to lower energies 
in fact requires additional features. For example 
a time resolved spectral study of GRB 041006 us- 
ing HETE-2 data covering the energy range of 
2 keV - 400 keV, indicated the presence of soft 
components, each of them having distinctive time 
evolution (Shirasaki et al. 2008). We emphasize 
again here that the use of different instruments can 
help in pinning down the systematics, identifying 
additional spectral components and developing an 
empirical spectral formalism. The simultaneous 
data from Swift /BAT and Fermi /GBM have the 
potential to explore the existence of additional 
components because the highly sensitive BAT has 
an energy range ('^15 keV - 200 keV) which is a 
subset of the GBM energy range. Our joint anal- 
ysis of BAT and GBM data has shown that in the 
20 - 40 keV region there is a distinct discrepancy. 
This could be due to a) systematic errors in each 
of the detectors, b) presence of additional compo- 
nents or c) an interplay of both these. Since the 
spectral fitting method used in XSPEC is forward 
in nature (that is fold the assumed spectrum and 
compare with the data - see Arnaud 1996), the 
derived spectral parameters will have strong bias 
of the bandwidth of the detectors, particularly for 



such low resolution spectral data. Hence a joint 
fitting of data with good overlap will result in a 
better estimate of the parameters. 

In conclusion, we have developed a method to 
test empirical descriptions of the energy spectra 
and their evolution with time in GRBs and to 
derive the underlying physical parameters. This 
is the first attempt in this different approach to 
study various parameters and predict quantities 
like spectral delay and pulse width of individual 
pulses in a GRB. 
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